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Abstract:  A  complete  characterization  is  given  of  closed  shift-invariant  subspaces  ofli^ffltrf) 
‘whldTprovide  a  specified  approximation  order.  When  such  a  space  is  principal'(i.e7,~generated  by  a 
single  function),  then  this  characterization  is  in  terms  of  the  Fourier  transform  of  the  generator.  As 
a  special  case,  we  obtain  the  classical  Strang-Fix  conditions,  but  without  requiring  the  generating 
function  to  decay  at  infinity.  The  approximation  order  of  a  general  closed  shift-invariant  space  is 
shown  to  be  already  realized  by  a  specifiable  principal  subspace.  r 
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1.  Introduction 


We  are  interested  in  the  approximation  properties  of  closed  shift-invariant  subspaces  of  the 
space  L2{WLd)  of  all  complex-valued  square-integrable  functions  on  ntd.  We  say  that  a  space  S  of 
functions  defined  on  is  shift-invariant  if,  for  each  /  €  5,  the  space  S  also  contains  the  shifts 
/(•  +  a),  a  €  TLA.  In  other  words,  <5  contains  all  the  integer  translates  of  /  if  it  contains  /.  A 
particularly  simple  example  of  such  a  space  is  provided  by  the  space 

M*) 

of  all  finite  linear  combinations  of  shifts  of  <j>.  We  call  its  /^(IR^-closure  the  principal  shift- 
invariant  space  generated  by  <f>  and  denote  it  by 

S(4>). 

Of  course,  a  general  closed  shift-invariant  subspace  of  £2(IR-d)  need  not  be  principal;  it  need  not 
even  be  generated  by  the  shifts  of  finitely  many  functions. 

Shift-invariant  spaces  are  important  in  a  number  of  areas  of  analysis.  Many  spaces,  encountered 
in  approximation  theory  and  in  finite  element  analysis,  are  generated  by  the  shifts  of  a  finite  number 
of  functions  <f>  on  IRd.  Shift-invariant  spaces  also  play  a  key  role  in  the  construction  of  wavelets.  In 
each  of  these  applications,  one  is  interested  in  how  well  a  general  function  /  can  be  approximated 
by  the  elements  of  the  scaled  spaces  Sh,  h  >  0,  of  S  defined  by 

Sh  :=  {s{-/h) :  s€S}. 

We  postpone  discussion  of  the  literature  concerning  the  approximation  by  the  scaled  spaces  Sh 
until  we  have  introduced  some  additional  terminology  and  stated  our  main  results. 

Associated  to  any  closed  subspace  S  of  ^(Hv  )  and  /  €  //2(Hld),  we  have  the  approximation 
error: 

(1.1)  £(/,£)  :=min{||/-s||:  s  €  S} 

where  here  and  throughout  this  paper  unsubscripted  norms  are  the  /^(DT^-norm.  We  are  interested 
in  describing  the  properties  of  S  which  govern  the  possible  decay  rates  of  E(f,Sh),  as  h  — ►  0,  (for  / 
in  j&2(IRd)  °r  in  one  of  the  smooth  subspaces  of  X2(n^-d))-  For  example,  we  shall  characterize  when 
the  scaled  subspaces  Sh,  h  >  0,  are  dense  in  the  sense  that  lim*-.o  E(f,Sh)  =  0.  More  generally, 
we  shall  characterize  when  the  spaces  Sh  approximate  suitably  smooth  functions  to  order  0(hk) 
as  h  — *  0. 

Our  definitions  of  approximation  orders  are  in  terms  of  the  potential  space  ^(IR'1),  k  >  0, 
defined  by 

wffltt')  :=  {/  €  i,(R') :  ll/ll, vJ(n.j  :=  (2ir)-'/J||(l  +  |  •  |)*/||  <  oo}. 

(Here  and  later,  we  use  |x|  :=  (x\  +  +  x^)1/2  t0  denote  the  Euclidean  norm  of  a  point  x  = 

(xi,...,Xd)  in  fit*1.)  When  A:  is  a  positive  integer,  these  are  the  usual  Sobolev  spaces.  We  say  that 
S  provides  approximation  order  k  if,  for  every  /  €  W2fc(lR<1), 

(1-2)  EU,S“)  <  c5/.‘||/||iv;(n',- 
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We  shall  characterize  the  closed  shift-invariant  spaces  S  which  have  this  property. 

A  variant  of  this  problem  is  to  characterize  when,  for  a  given  k  >  0,  we  have  for  each  /  € 
W£(m.d)  (in  addition  to  (1.2)), 

(1.3)  E(f,Sh)  =  o(hk),  0. 

When  Jb  =  0,  this  is  just  the  density  problem.  For  this  reason,  we  say  that  S  provides  density 
order  k  whenever  (1.3)  holds. 

Our  characterizations  of  density,  approximation  order,  and  density  order  are  in  terms  of  Fourier 

j  ^ 

transforms.  If  /  €  ii(IR  ),  its  Fourier  transform  /  is  defined  by 

f(y)  :=  /  f{x)e-ix-y  dx. 

Jr.* 

We  assume  that  the  reader  is  familiar  with  the  usual  properties  of  Fourier  transforms.  In  particular, 
the  Fourier  transform  can  be  uniquely  extended  to  l2(ntd)  and  more  generally  to  the  space  of 
tempered  distributions  on  IR. 

Many  authors  have  shown  (under  various  restrictive  conditions  on  <j>)  that  the  approximation 
properties  of  a  principal  shift-invariant  space  S(<t>)  are  related  to  the  order  of  the  zeros  of  the  Fourier 
transform  of  <f>  at  the  integer  multiples  of  2ir.  It  is  therefore  not  surprising  that  our  characterizations 
of  approximation  order  for  general  principal  shift-invariant  spaces  involve  the  behavior  near  zero 
of  the  2jr-periodization  of  \<p\2,  i.e.,  the  L2(TT ^-function 

(i-4)  i?.a »  e 

/3€2ir22<' 

This  function  enters  our  considerations  as  part  of  the  function  6  Loo(C),  defined  on  the  cube 

C  :=  [— 7T . .  7r]d 


by 

a  2 

(1.5)  A,  :=  (1  -  on  C. 

Here  (and  below  without  further  comment),  we  make  use  of  the  obvious  identification  between  the 
space  L2(T£d)  of  functions  on  the  d-dimensional  torus  'H’d  and  the  space  L2(C)  of  functions  on  the 
fundamental  domain  C.  We  shall  prove 

Theorem  1.6.  The  principal  shift-invariant  subspace  S(<f>)  of  ^(IR^)  provides  approximation 
order  k  >  0  if  and  only  if  |  ■  |-*A$  is  in  L 00(C). 

The  analogue  of  this  result  for  density  orders  is 

Theorem  1.7.  The  principal  shift-invariant  subspace  S(<f>)  of  T2(lRd)  provides  density  order  k  >  0 
if  and  only  if  |  •  |-feA0  is  in  L 00(C)  and 

(1.8)  Jim  A-"  /  |»J-**[A*C»)|*  =  0. 

"-*0  JhC 

Of  course,  in  the  case  k  =  0,  (1.8)  characterizes  when  we  have  density. 

It  is  rather  remarkable  that  these  conditions  also  characterize  approximation  and  density  orders 
for  arbitrary  closed  shift-invariant  subspaces  of  Z2(IRd)'  Namely,  we  shall  prove: 
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Theorem  1.9.  A  closed  shift-invariant  subspace  S  of  Z>2(IR,d)  provides  approximation  order  k  >  0 
if  and  only  if  it  contains  a  function  <j>  for  which  |  •  |-feA^  is  in  L^C).  The  space  S  provides  density 
order  k  >  0  if  and  only  if  it  contains  a  function  4>  for  which  |  •  €  Zoo(C)  and  (1.8)  holds. 

We  prove  the  last  theorem  by  showing  in  §3  that  the  case  of  approximation  by  arbitrary  closed 
shift-invariant  subspaces  of  Z^IR1*)  can  be  reduced  to  the  case  of  principal  shift-invariant  spaces. 

In  the  case  of  principal  shift-invariant  spaces,  our  method  of  proof  is  based  on  two  results 
which  we  feel  will  have  other  important  applications.  The  first  is  an  explicit  formula  for  the  best 
Z2(IR<1)-approximation  from  S(<f>).  The  second  is  the  following  characterization 

(1.10)  S(<f>)  =  {t4>  €  Z2(IRrf) :  t  is  27r-periodic} 

of  the  space  S(<f>)  in  terms  of  its  Fourier  transform.  Here  and  later,  for  a  set  F  of  functions,  we 

A  A. 

denote  by  F  :=  {/  :  /  6  F}  the  set  of  its  Fourier  transforms. 

It  turns  out  that  our  analysis  applies  equally  well  to  the  more  general  situation  where  the 
h- refinement  of  the  space  S  is  obtained  by  means  other  than  scaling.  Such  cases  are  known  and  are 
of  interest  in  both  spline  theory  (e.g.,  exponential  box  splines,  cf.  [DR])  and  radial  basis  function 
theory  (cf.  the  detailed  discussion  in  [BR2]).  In  the  non-scaling  case,  we  employ  a  family  {Sh}h  of 
shift-invariant  spaces,  and  consider  the  rates  of  decay  of  E(f,Sfc)  as  a  function  of  h.  The  notions 
of  “approximation  order  k"  or  “density  order  fc”  for  the  sequence  {£»,}/,  are  obtained  by  replacing 
each  E(f,Sh)  in  the  above  definitions  by  E(f,Sfc). 

We  close  this  section  with  a  brief  discussion  of  the  connections  between  the  results  of  this 
paper  and  results  in  the  literature.  Schoenberg,  in  his  seminal  paper  [S],  was  the  first  to  recognize 
the  importance  of  the  Fourier  transform  for  describing  approximation  properties  of  principal  shift- 
invariant  spaces.  For  the  case  d  =  1,  and  with  <j>  a  piecewise  continuous  function  with  exponential 
decay  at  infinity,  Schoenberg  showed  that  all  algebraic  polynomials  of  degree  <  k  can  be  written 
in  the  form  YlatTi'*  <K’  ~  a)c(a)  in  case 

(1.11)  <£(0)  /  0  and  Dy$  =  0  on  2x7Zd\0  for  all  I7I  <  k 
holds  (with  d  =  1). 

Strang  and  Fix  [SF]  have  treated  the  approximation  properties  of  the  space 

of  all  linear  combinations  <f>( •  -  a)c(a)  (finite  or  not)  of  the  integer  shifts  of  a  compactly 

supported  function  </>.  There  is  no  problem  of  convergence  of  such  sums  since,  for  any  point  x  €  IR**, 
at  most  finitely  many  terms  of  the  sum  are  nonzero  at  x.  Strang  and  Fix  necessarily  restricted 
attention  to  the  subspace 

S2{<f>)  :=  «S*(4>)  n  Z2(IRd). 

While  this  space  is,  in  general,  not  closed  in  Z2(lRd),  one  can  show  (see  Theorem  2.16  below) 
that  its  Z2(IR<<)-closure  is  <S(d>).  Strang  and  Fix  proved  that  52(<£)  provides  approximation  order  k 
whenever  (1.11)  holds. 
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To  compare  this  result  with  Theorem  1.6  above,  note  that,  for  a  compactly  supported  <f> ,  [4>,4>] 
is  a  trigonometric  polynomial,  since  then 

(1.12)  [<f>,<j>\  =  ^2  a(a)ea,  with  a(a)  :=  j  <f>(x  -  a)<f>{x)dx. 

JlRd 

Here  and  later,  we  use  the  abbreviation 

ea(y)  :=  ete*». 

If  (1.11)  holds,  then  [<£,  <f>]  does  not  vanish  at  the  origin  and  of  (1.5)  has  a  zero  of  multiplicity  k 
there.  Thus,  |  - 1~ is  in  Loo(C)  (as  we  know  it  must  be).  However,  there  are  two  important  points 
to  bear  in  mind  concerning  our  Theorem  1.6  and  the  Strang-Fix  result.  First  of  all,  our  theorem 
does  not  require  that  <f>  be  compactly  supported,  nor  even  that  it  decay  at  infinity.  Secondly,  it 
applies  even  when  <f>  vanishes  at  the  origin,  a  case  of  practical  importance  yet  not  accessible  to 
earlier  approaches. 

Actually,  Strang  and  Fix  proved  more  than  we  have  just  stated  since  they  showed  that  the 
approximation  order  0(hk)  to  a  given  /  €  by  the  elements  of  can  be  achieved 

with  a  control  on  the  coefficients  of  the  approximants  €  S2(<j>)k.  Namely,  if  the  approximants 
are  represented  with  respect  to  the  X2-normalized  functions  <f>(a,h,  x)  :=  h~dI2<f>(x/h  —  a)  by 
Sh  =  £<*€2*  c/>(«)<K«i  ')> then 

(1-13)  ll(c*(a))IUa(*-)  <  cf- 

The  introduction  of  such  controlled  approximation  is  important,  since  Strang  and  Fix  show 
that,  conversely,  if  5a(«f>)  provides  controlled  approximation  order  k  then  (1.11)  holds.  In  other 
words,  for  compactly  supported  <f>,  S2(<f>)  provides  controlled  approximation  order  k  if  and  only  if 
(1.11)  holds.  Since  it  can  be  easily  seen  that  our  condition  in  Theorem  1.6  is  weaker  than  (1.11) 
(even  for  compactly  supported  <f>),  it  follows  that  there  are  cases  when  the  achievable  approximation 
order  cannot  be  obtained  in  a  controlled  manner.  In  this  connection,  it  is  worthwhile  to  point  out 
that  positive  controlled  approximation  order  forces  <£(0)  ^  0. 

There  is  a  rich  literature  of  clarifications  and  extensions  of  the  Strang-Fix  result,  including 
extensions  to  noncompactly  supported  <f>  ([BH2],  [J2],  [DM2],  [BJ],  [Bl],  [R],  [CL],  [JL],  [HL] 
[BR2]).  In  addition,  there  are  many  papers  studying  the  approximation  order  of  specific  principal 
(and  other)  shift-invariant  spaces,  some  of  them  ([Bui, 2],  [BD],  [BuD],  [Bill],  [BR1],  [DJLR], 
[DM1],  [DR],  [Ja],  [J  1],  [L],  [LJ],  [M],  [MN1,2],  [Ra],  [RS])  are  included  in  the  references;  see 
also  the  surveys  [B2],  [C],  [P]  and  the  references  therein.  By  making  assumptions  on  <j>  weaker 
than  those  used  in  any  of  the  above  references,  we  can  translate  our  conditions  on  A^  into  simple 
conditions  on  4>.  For  example,  we  show  in  §5  the  following: 

Theorem  1.14.  Assume  that  <t>  is  bounded  on  some  neighborhood  of  the  origin.  If  provides 
approximation  order  k,  then  <f>  has  a  zero  of  order  k  at  every  (3  €  2ir7Ld\0.  In  particular,  ZP  <£(/?)  =  0 
for  all  |7|  <  k  in  case  <f>  is  k  times  differentiable  (in  the  classical  sense)  at  such  (3. 

Note  that  the  boundedness  of  <f>  required  here  holds,  for  example,  if  <j>  is  continuous  at  0.  In 
particular,  it  holds  for  every  4>  €  Xi(Dld). 

We  also  show  in  §5  the  following  converse: 
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Theorem  1.15.  Assume  that  1/0  is  bounded  on  some  neighborhood  of  the  origin  and  that,  for 
some  p  >  k  +  d/2,  all  derivatives  of  </>  of  order  <  p  are  in  12(A),  with  A  :=  Be  +  (27rZd\0)  for  some 
open  ball  Be  centered  at  the  origin.  If  Dy4>(P)  =  0  for  all  I7I  <  k  and  all  (3  €  2ir7Zd\0,  then  S(<f>) 
provides  approximation  order  k. 

For  most  of  the  examples  of  a  non-compactly  supported  </>  in  the  literature  (e.g.,  radial  basis 
functions,  see  [P]),  <f>  is  very  smooth  on  IRd\0,  but  has  a  singularity  at  the  origin.  On  the  other 
hand,  the  present  standard  approach  to  the  derivation  of  approximation  orders  (viz.,  the  polynomial 
reproduction  argument)  requires  <f>  to  decay  at  00  (at  least)  like  0( |  •  |"(fc+d)),  hence  requires  <f>  to 
be  globally  smooth.  To  circumvent  this  obstacle,  one  usually  seeks  a  function  0  €  «So(0)  (or  in 
some  superspace  of  So(<t>))  whose  Fourier  transform  0  is  smoother  than  0,  since  this  implies  a  more 
favorable  decay  of  0  at  00.  This  ‘localization’  process  constitutes  the  main  effort  in  establishing 
approximation  orders  for  a  non-compactly  supported  4>.  Our  theorem,  though,  does  not  require  <f> 
to  decay  at  00  at  any  particular  rate,  thus  obviating  the  search  for  such  0.  Results  (weaker  than 
the  above  theorem)  about  Xoo(IRd)-approximation  orders,  that  apply  to  functions  which  decay 
only  mildly  at  00,  were  derived  in  [BR2].  The  approach  there  exploits  the  fact  that  the  exponential 
functions  e$,  0  6  IRd,  are  in  the  space  in  which  approximation  takes  place.  In  contrast,  the  approach 
here  makes  use  of  the  simple  and  explicit  formula  for  the  orthogonal  projection  onto  <S(0). 

2.  The  orthogonal  projector  onto  5(0) 

In  this  section,  we  derive  two  important  facts  about  the  principal  shift-invariant  space  5(0) 
which  will  be  the  basis  of  much  of  the  analysis  that  follows.  The  first  is  a  simple  formula  (given  in 
Theorem  2.9)  for  the  (Fourier  transform  of  the)  best  /^-approximation  from  5(0).  The  second  is 
the  description 


(2.1) 


5(0)  =  {r0  €■  l2(IRd)  :  T  »s  27r-periodic) 


of  S(<j>)  in  terms  of  Fourier  transforms  mentioned  in  the  introduction. 

The  yet  to  be  proven  (2.1)  suggests  that  the  calculation  of  integrals  and  inner  products  in¬ 
volving  functions  from  S(<t>)  should  be  taken  over  the  torus  Trd.  This  can  be  accomplished  by 
periodization.  If  g  €  Li(IRd),  then 


(2.2) 


/  9=  £  /  9=  f  g\ 
J*d  P^Jc^  Jc 


with 

g°~  £  j(-  +  /3) 

0€2w  2Z* 

the  (2x)-periodization  of  g.  Here,  the  sum  is  to  be  taken  in  the  sense  of  £i('ITd)-convergence, 
which  makes  sense  since,  by  assumption,  g  €  Ii(IRd).  In  particular,  g°  6  Ii('ITd). 

Similarly,  we  have 


(2.3) 


/  9o9\  =  /  [50,01] 
J  n*  Jc 
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for  the  inner  product  of  two  functions  go,gi  €  £2(1R*),  with 


(2.4)  [ffo>0i]  :=  (0o<7i)°  =  0°('  +  0)9i(-  +  /?)• 

/}€2irEJ 

Note  that  lgo,gt]  is  in  Xi('IT<1)  since  goffi  €  Li(JRd).  Also,  by  the  Cauchy- Schwarz  inequality, 

(2.5)  |[<7o,<7i]|2  <  [<?o,<7o][<7i,fli], 


and  the  right  side  of  (2.5)  is  finite  a.e.  We  will  most  often  use  (2.3)  in  the  form 


§  ^  i 

which  is  valid  for  arbitrary  f,<f>  €  X2( IR  )  and  arbitrary  27r-periodic  r  for  which  r/  €  L2(IR  ).  We 
note  that  (2.6)  implies  the  estimate 

(2-7)  Ik^IlLatBV1)  ^  IMl£2{Trd)ll[0>0]llLoo(’ir<<) 

of  use  when  [<j>,  4>\  is  bounded,  e.g.,  when  4>  is  compactly  supported. 

After  these  brief  remarks,  let  us  consider  the  problem  of  finding  a  formula  for  the  projection 
of  L2(IRd)  onto  S(<f>).  Let  P  :=  P#  denote  the  orthogonal  projector  onto  S(<j>).  Then  Pf  is  the 
unique  best  L^IR^-approximation  to  /  from  S(<t>),  and  is  characterized  by  the  fact  that  it  lies 
in  S(<f>)  while  its  difference  from  /  is  orthogonal  to  Since  the  Fourier  transform  preserves 

orthogonality,  it  follows  (for  example  from  the  uniqueness  of  best  approximation  in  L2(1R“))  that 
the  orthogonal  projector  P  onto  S(<f>)  satisfies  Pf  =  Pf. 

We  consider  first  what  it  means  for  a  function  /  to  be  orthogonal  to  S(<f>).  Since  finite  linear 
combinations  of  the  (integer)  shifts  <f>{ •  +  a)  of  0  are  dense  in  S{<}>),  f  €  L2(flld)  is  orthogonal  to 
S(<{>)  iff  /  is  orthogonal  to  e_Q<£  for  every  a  6  Z5d,  i.e.  (with  (2.6)),  iff 

0  =  /  /  eQ<j>  =  j  [J,4>]ea  for  all  a  €  5Zrf. 

J  it*  Jc 


This  proves 

Lemma  2.8.  The  orthogonal  complement  S(^)1  of  S(<f>)  in  L2(Dl<i)  consists  of  exactly  those 
f  €  L2( IRd)  for  which  [fjj  =  0. 

From  Lemma  2.8,  we  can  easily  determine  Pf.  Suppose,  as  is  suggested  by  (2.1),  that  Pf  =  r<p, 
with  r  some  2^-periodic  function.  Then,  from  Lemma  2.8, 

[7»^]  =  [Pfi<f>]  =  [t4>,4>]  —  t[4>,4>\. 


This  motivates  the  following: 
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_  1  — •  ™  ^  ^ 

Theorem  2.9.  For  each  f  G  Z/2(IB  ),  P<pf  =  Tj<f>,  with  the  2ir -periodic  function  tj  defined  by 
(2.10)  iy:=| 


L/\<£]/[<M]  on 
0  otherwise, 


and  flj,  defined  up  to  a  null-set  by 

ft*  :=  supp$,$]  :=  {u  G  TTd  :  $,£](a/)  #  0}. 

Proof.  It  is  enough  to  show  that  Pf  =  Tj<f>  for  each  /  G  I-2(ni  ).  We  first  want  to  see 
that  Tj<f>  is  in  L2{fR.d).  By  (2.5),  \Tf\2[<f>,<f>]  <  [/,/}.  With  this,  two  applications  of  (2.6)  give 

(2.11)  /  Ml ’=  /  MlM< /l/.7)=  /  l/lJ- 

./ii'  Jc  Jc  Jr* 

Consequently,  Tj<f>  G  f/2(IR)  and  moreover  the  linear  map 

Q  :  ^(ffi-11)  -  L2(m.d)  :f~Tf4> 

is  well-defined  and  norm- reducing  on  X2  (IR.^).  We  next  prove  that  Q  =  P. 

If  /  G  S(<j>)  =  (S(<f>)Ly,  then  Lemma  2.8  gives  that  77  =  0,  hence  Qf  =  0.  Thus  Q  =  P  on 
S(<f>)  -.  On  the  other  hand,  on  SI#  =  supp[<£,$, 

=  l ea4>AV\4>,4>]  -  «o,  for  all  a  G  TLd . 

Since  ^  =  0  on  the  complement  of  +  27rZ5  ,  this  implies  that  Q  maps  the  Fourier  transform 
of  each  integer  shift  of  <f>  to  itself.  Since  Q  is  linear  and  bounded,  and  coincides  with  P  on  a 
fundamental  set  for  S(<j>),  we  have  Q  =  P  on  S(<f>).  By  linearity,  Q  =  P  on  all  of  L2(IR-‘<)-  ♦ 

Remark.  With  the  convention  (which  we  use  throughout  this  paper)  that  0  times  any  ex¬ 
tended  number  is  0,  we  are  entitled  to  write 


(2.12) 


TJ  =  [/>$/[<£»£]  and 


p~~f  _  [/>$-? 


Note  that  (2.11)  supplies  the  following  lemma. 

Lemma  2.13.  If<f>,f  G  L2(JR.d),  then  T/<j>  G  ^2(IILd),  and  \\tj4>\\  <  ||/||. 

As  a  consequence,  we  obtain  the  characterization  (2.1)  of  the  space  S{<f>)  in  terms  of  its  Fourier 
transform. 

Theorem  2.14.  A  function  f  is  in  S(<f>)  if  and  only  if  f  =  r<f>  for  some  2z -periodic  r  with 
T<p  G  J^OR-  )•  fo  particular,  r<p  G  S(4>)  for  every  bounded  r. 

Proof.  If  /  G  S(<f>),  then  Pf  =  f.  Hence,  by  Theorem  2.9,  /  =  with  77  the  2ff-periodic 
function  [/,$)/$,$),  and  Tj<f>  G  £.2(IR‘f )  because  of  Lemma  2.13. 
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a  ^  j 

Conversely,  if  r  is  defined  on  TT  ,  and  r<f>  6  Z2(IR),  then  the  inverse  transform  /  of  r<f>  is  also 
in  L2(JSLd)  and  satisfies  r/  =  [r<£,  <£)/[<£,  7]  =  r  on  fl*  =  supp[7,#  Since  4>  vanishes  off  Q^  +  2?rZZrf, 
this  implies  with  Theorem  2.9  that 

Pf  =  tj4>  =  T<j>  =  f. 

Consequently,  Pf  =  /  and  hence  /  €  <S(<£).  Finally,  if  r  is  bounded,  then  r<f>  €  Z2(IR)  since 
4>  €  Z2(IR<1).  4 

Remark  2.15.  Asher  Ben-Artzi  has  pointed  out  to  us  that  Theorem  2.14  could  have  been 
derived  from  general  results  (cf.  Theorem  8  of  [H;p.59])  concerning  closed  subspaces  of  ^(TT,^) 
which  are  invariant  under  multiplication  by  exponentials.  Furthermore,  the  lemma  of  [H;p.58] 
shows  that  Theorem  2.14  implies  Theorem  2.9. 

Remark.  The  representation  r<f>  for  /  6  S(<j))  is  in  general  not  unique.  If  t04>  =  rx<£,  we  can 
only  conclude  that  r0  =  T\  a.e.  on  £1$.  However,  if  the  shifts  of  <j>  are  an  orthonormal  basis  or,  more 
generally,  an  Z2(IRr)-stable  basis,  then,  as  is  well  known,  [4>,  0]  and  its  reciprocal  are  both  in 
and  not  only  is  the  representation  unique  but  the  function  r  is  in  Z2('ITd).  It  is  interesting  to  note 
further  that  we  have  a  unique  representation  even  when  the  shifts  of  4>  are  not  an  Z2(lR<1)-stable 
basis  provided  £1$  differs  from  H’<<  only  by  a  null-set. 

The  following  consequence  of  Theorem  2.14  is  of  importance  when  comparing  our  results  with 
related  results  in  the  literature. 

Theorem  2.16.  If  4>  €  I>2(Htrf)  has  compact  support,  then  S{4>)  is  the  L2{JRd)-closure  of  52(^>)  = 
S*(*)nz2(lRd). 

Proof.  Since  S(<f>)  is  the  Z2(IR‘1)-closure  of  So{<f>)  and  So(<f> )  is  contained  in  S2(<t>)  (since 
</>  6  Z2(IR<<)),  we  only  have  to  prove  that 

(2.17)  S2(j)  C  S(<t>). 

We  now  prove  this  by  showing  that  P#f  =  /  for  every  /  G  S2(0),  i.e.,  with  (2.12),  that 

(2.18)  /  = 

Since  <f>  has  compact  support,  [4>,<t>]  is  a  trigonometric  polynomial  (cf.  (1.12)),  hence  (2.18)  is 
equivalent  to  the  equation 

(2.19)  [$,<!>]  f  =  [/,0]0  a.e., 

and  it  is  this  equation  we  now  verify  for  any  /  in  Z2(IR‘i)  of  the  form  'ShpeTL4  00  ~  0)c(0)- 
We  do  this  by  showing  that  both  sides  of  (2.19)  are  the  Fourier  transform  of  the  function 

Y  /(  +or)a(o), 
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with  a(a)  =  JRd  <f>( — a)<f>  the  (Fourier)  coefficients  of  the  trigonometric  polynomial  [<j>,  <f>],  see  (1.12). 
This  is  immediate  for  the  left  side  of  (2.19)  since  (Ya$7zd  /(•  +  <*)a(a))T  =  (S a<z"Zd  a(a)ea)f  for 
any  /  €  Z2(IRd)  and  any  finite  sequence  (a(a)),  and  [4>,4>]  is  indeed  a  finite  sum  of  exponentials 

A  A 

since  <f>  is  compactly  supported.  As  to  the  right  side  of  (2.19),  [/,  <f>\  is  a  2x-periodic  Z,2-function 

^  j 

(since  <f>  is  compactly  supported,  thus  <f>  is  bounded),  hence  the  ^(Tr )-limit  of  its  Fourier  series 
E7g nd  K"i)eT>  with  b  given  by 

b(n)  ■=  (2n)~d  f  [f,$\e.y  =  j  /^(-  +  7)=  £  [  <t>{-  -  PW)<K-  +  7)  =  J]  #(7  +  /J). 
*c  *'1R,<I  0&rai  *'R<<  p£iLd 

By  (2.7),  [J,4>)4>  is  the  I,2(IRd)-liinit  of  Y^en,d  &(7)e7^»  hence  ([/,$]0)v  is  the  I2(lRd)-limit  of 
'H-fe2fr>  ^('  +  7)^(7)*  Since  this  last  sum  also  converges  uniformly  on  compact  sets,  these  two  limits 
must  be  the  same.  This  implies  that  the  right  side  of  (2.19)  is  the  Fourier  transform  of 

^(•  +  Tr)  c(/?M7  +  £)=  J]  <K- +  *  - /3)c(/3)a(a)  =  /(•  +  a)«(a), 

•>€  2Z*1  tfgZ''  aeZZ11  aez*1 

with  the  rearrangement  of  the  sums  justified  by  the  fact  that  all  sums  are  finite.  4 

We  now  turn  to  our  main  objective,  viz.  the  error  of  the  best  approximation.  If  /  is  supported 
in  one  the  cubes  f3  +  C,  (3  €  2ir7ld,  this  error  takes  a  very  simple  form: 

Theorem  2.20.  Let  <f>  €  I.2(IRd).  If  /  €  Z,2(lRd)  and  supp  /  C  (3  +  C  for  some  (3  €  2tt7Ld,  then 

(2.21)  =  (2^)-'£l(7,2f(5)):!  =  {2n-)-<'  /  l/|2(l-#W). 

j  R*  [<M] 

Proof.  Since  supp  /  C  C  +  (3  for  some  /3  €  27r5Zd,  we  have  [/,<£]  =  /(•  +  /3)<t>( ■  +  (3)  on  C. 
Therefore,  with  (2.6), 

imii2=  /  I  \f\2\t\2/ii*)- 

Jc  Jlid 

By  Theorem  2.9,  this  shows  that 

=  (2*)"'  /  l/|2|5l2/(?,?l, 

Jnd 

and  this  finishes  the  proof  since  ||/  -  P^/j|2  =  ||/||2  -  ||P^/]|2.  4 

3.  The  reduction  to  the  principal  case 

The  explicit  and  simple  expression,  derived  in  the  previous  section,  for  the  orthogonal  pro¬ 
jector  onto  a  principal  shift-invariant  space  will  also  prove  to  be  very  useful  in  the  discussion  of 
approximation  from  a  general  shift-invariant  space.  For,  remarkably,  the  approximation  power  of 
a  general  shift-invariant  space,  however  large,  is  already  contained  in  a  single  (suitably  chosen) 
principal  shift-invariant  subspace  of  it.  The  next  proposition  provides  the  algebraic  background  for 
this  fact.  We  use  repeatedly  the  simple  observation  that  the  best  approximation  Pf  to  /  from  S 
is  also  the  best  approximation  Ppjf  to  /  from  S(Pf ),  i.e., 

Ppjf  =  Pf. 
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Proposition  3.1.  Let  P  be  the  orthogonal  projector  onto  the  closed  shift-invariant  subspace  S  of 
i2(IRa)  and  denote  by  P  the  corresponding  orthogonal  projector  onto  S.  Then  P(r/)  =  rP(f)  for 
any  f  G  /^(IR*1)  and  any  2r -periodic  r  for  which  rf  G  X2(IR‘1). 

Proof.  If  S  is  principal,  then  the  conclusion  follows  directly  from  (2.12).  For  the  general 
case,  the  assumptions  on  r  and  /  imply  with  Theorem  2.14  that  rf  G  S(f).  Since  S(f)  is,  by 
definition,  the  F2(lR“)-closure  of  5o(/),  and  So(f)  =  {rn/  :  rn  a  trig.polynomial},  it  follows  that  rf 
is  the  P2(Ht  )- limit  of  rn/  for  some  sequence  (rn)  of  trigonometric  polynomials.  The  shift-invariance 
of  S  and  the  uniqueness  of  the  best  /^-approximation  imply  that  P(/(--f  a))  =  (P/)(--fa)  for  every 
/  G  i2(IR<<)  and  every  a  G  7Ld.  Hence,  taking  finite  linear  combinations  of  Fourier  transforms, 
P(rn/)  =  rnP/,  and  so,  by  the  continuity  of  P, 

P(rf)  =  lim  P(r„/)  =  lim  rnPf. 

n— *•  oo  n— ►  oo 

Each  r„P/  is  in  the  closed  space  S(Pf),  therefore  also  P(r/)  lies  in  S(Pf).  Thus,  projecting  rf 
onto  S  is  the  same  as  projecting  it  onto  the  subspace  S(Pf)  of  S.  Since  we  already  know  that 
P^(r/)  =  rP^f  for  any  <j>,f  G  Z2(IR<i),  this  means  that  we  obtain 


Hr?)  =  Ppf(rf)  =  rPPI(f)  =  rP/, 


the  last  equality  since  PPjf=  Pf.  4 

Corollary  3.2.  If  P  is  the  orthogonal  projector  onto  some  shift-invariant  subspace  of  £2(IRd)  and 
g  G  X2(IR<1),  then 

PPg  =  PpgPg. 

Proof.  If  /  G  Z/2(md),  then  Pgf  =  rg  for  some  27r-periodic  r  and  therefore  by  Proposition 
3.1,  P{Pgf)  =  rPg.  On  the  other  hand,  Pp9(P9f)  =  Ppg{Tg)  =  rPpag  =  rPg.  4 

Theorem  3.3.  For  any  closed  shift-invariant  subspace  S  of  F2(lRrf)  and  any  f,g  G  Z-2(IRd), 

(3.4)  E(f,S)<  E(f,S(Pg))<  E(f,S)+2E(f,S(g))1 

with  P  =  Ps  the  orthogonal  projector  onto  S. 

Proof.  Only  the  second  inequality  needs  proof.  By  Corollary  3.2, 

f-Pp9f  =  f-Ff  +  Pf-PP9f  +  Pp9P9f-Pp9f , 

and  therefore 

(3-5)  11/  -  PP9f\\  <  ||/  -  P/ll  +  ||/  -  P,/||  +  ||P5/  -  /||. 
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This  theorem  shows  that  the  approximation  order  of  the  particular  principal  subspace  S(Pg ) 
of  S  is  the  same  as  that  of  all  of  S ,  provided  that  the  approximation  o'der  of  the  principal  space 
S(g)  is  at  least  as  good  as  that  of  S.  This  suggests  the  use  of  a  special  function  g *  for  which  S(g*) 
has  arbitrarily  high  approximation  order.  We  can  take  g *  to  be  the  inverse  Fourier  transform  of 
the  characteristic  function  of  the  cube  C  =  [— ir .  .ir]4,  i.e., 

5*  :=  (XC)V- 

We  note  that,  by  (2.12),  fP~f  =  [J,XC\I\XC<XC  1  Xc  =  Xc7-  Hence, 

(3.6)  E(f,S(g-))  =  (2jr)-‘*/2[|(l  -  xc)7l 

This  allows  us  to  show  easily  that  the  space  S(g’)  provides  approximation  and  density  order  k  for 
all  k  >  0.  For  this,  we  follow  the  example  of  [BR2]  and  consider,  equivalently,  the  approximation 
of  the  scaled  function 

h  ==  /(/»•) 

from  the  fixed  space  S  instead  of  the  approximation  of  the  function  /  from  the  scaled  space  Sh. 
For, 

(3.7)  E(f,Sh)  =  hd'2E(fH,S), 


as  is  easily  established  by  a  change  of  variables. 
Lemma  3.8.  Let  f  6  W^nt**),  k  >  0,  h  >  0.  Then 


E(/,S(g-)h)  <  </(/.)** l/llwsi*). 

with  the  (nonnegative)  function  Cf  defined  by 


(3.9)  ej(hf 

hence  €/(h)  <  1,  and  C/(0+)  =  0. 


/(nd\c)/fc(*  +  I  ’  l)2fcl/l2 

/b-(i  +  H)**I/I*  ” 


Proof:  Since  /  €  W^IR*1),  the  function  v  :=  (1  +  |  •  |)*/  is  in  L2(®-<1)i  and  ||/||w*(r‘)  = 
(2jr)“rf/2||i/||.  Since  //,  =  h~df(-/h),  (3.7)  and  (3.6)  imply  that 


(2  *)iE(f,S{g-)k)2  =  (2*  )'*"£(/»,  SO,*))1  =  *'|(1  -  xc)hf  =  h‘  [  Mgtfdg 

Jtid\C 


-h'd  f  \f(y/h)\2dy=  h2k~d  I 
Jnd\c  Jt\ 


Hy/h)\7 


lnd\c 

<  h2k~d  f  W(y/h)\2  dy  =  h2k  [  \v\2  =  (2*)dh2kef(h)2\\f\\2wi (R-). 

Jnd\c  J{nd\c)/h  ’ 


R  d\C 
dy 


R d\c  +  M)2fc 

2  _  (OnrWfc2*,,  ./),\2||  ri|2 


We  note  for  later  reference  the  following  useful  result  established  during  the  proof  of  Lemma 


3.8. 
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Corollary  3.10.  For  each  f  €  W'''2*(Et<,)> 

*d/2||(l  -  XC)A||  <  V*)d,2cf(h)hk\\f\\wnn<‘), 

with  €/  given  by  (3.9). 

Now  let  S  be  an  arbitrary  closed  shift-invariant  subspace  of  and  let  <f> *  :=  Pg*  be  the 

best  i2(Eld)-approximation  to  g *  from  S.  Using  (3.7)  and  Lemma  3.8  in  (3.4),  we  obtain 

(3.11)  EU,Sh)  <  £(/,«(*•)*)  <  EU,Sk)  +  2£/(1.)A‘II/I|w.,r.,, 

with  €/(h)  given  by  (3.9).  This  means  that  S  provides  approximation  order  k  >  0  or  density  order 
fc  >  0  if  and  only  if  its  principal  shift-invariant  subspace  S(<f>*)  does.  More  than  that,  since  e/(h) 
does  not  depend  on  S ,  it  proves  the  following: 

Theorem  3.12.  The  sequence  {Sh}h  of  closed  shift-invariant  subspaces  of  L2(l Rd)  provides  ap¬ 
proximation  order  k  >  0  or  density  order  k  >  0  if  and  only  if  the  corresponding  sequence  (S(<f>^)}h 
of  principal  shift-invariant  subspaces  ( with  4>'h  :=  PsK(g*)  and  g *  =  \wc)  does. 

4.  Approximation  orders  and  density  orders 

In  this  section  we  give  a  complete  characterization  of  approximation  orders  and  density  orders 
from  the  sequence  {S^jh  of  shift-invariant  spaces.  In  view  of  Theorem  3.12,  we  need  only  to  consider 
the  special  case  when  each  Sh.  is  principal.  For  <f>  €  L2(Eld),  we  let  €  Loo(C')  be  defined  as  in 
the  introduction: 

a*:=(i'S),/j’  °"c- 

In  terms  of  this  A^,  (2.21)  gives  that 

(4.1)  E{f,Sm  =  (27r)-rf/2||/A0H  if  supp/  C  C. 

For  /  G  L2(Etrf)  with  /  not  just  supported  in  C,  we  estimate  E(f,S(<t>))  =  (2 x)~d/2E(f,S(4>)) 
with  the  aid  of  Corollary  3.10  and  the  simple  observation  that 

-  F(XC/J)\  <  11(1  -  Xc)/ll 

for  an  arbitrary  subspace  S  of  Z/2(Eld).  Indeed,  with  the  aid  of  (3.7),  this  estimate  implies  that 

!«(/,$*)  -  W2x)d/*E(XCTH,$)\  =  | h^EUk.S)  -  W2*)d'’E{xchJ)\ 

=  (h/2x)W\E(U,S)  -  E(Xch,S)\  <  (fc/2»)^||(l  -  XC)AI|. 

Therefore,  Corollary  3.10  establishes 

(4.2)  l^/.S1-)  -  (lt/2x)'"!£(xcA.S)l  <  </(W*ll/ll  «.*<«<)• 
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Theorem  4.3.  For  {4>h}h  C  L2(JR.d),  the  sequence  {S(4>h)}h  provides  approximation  order  k  if 
and  only  if 


{ 


■> 

(h  +  \-\)kih 


is  bounded  in  L0 0  ( c )• 


Remark.  Since  each  A^h  is  non-negative  and  bounded  above  by  1,  and  since  each  (h  +  |  •  |)fe  is 
bounded  below  by  hk,  it  is  clear  that  each  is  an  element  of  Loq(C).  So  it  is  the  uniform 

boundedness  of  as  h  — ►  0  that  characterizes  the  approximation  order  k. 

Proof.  In  view  of  (4.2),  {S(4>h)}h  provides  approximation  order  k  if  and  only  if  there  exists 
some  constant  c  such  that  for  every  /  6  W* (ntd)  and  every  h  >  0 

(4.4)  hinE(XcK,SM)  <  cA*||/|||v*(rJ). 

Since  Xcfh  is  supported  in  C ,  we  may  appeal  to  (4.1)  (i.e.,  to  Theorem  2.20)  to  conclude  that 


(4.5) 


hdE(Xch,SM?  =  hd  /  lAI’A^ 

J  c 


=  k-d  I  = 

JC 


For  /  e  W^nF*),  the  function  u  :=  (1  + 1  •  \)kf  is  in  /^(Dt"1)*  and  ||/||w‘(ir4)  =  (27r)  d/2||i/||.  With 
the  aid  of  i/,  the  last  expression  in  (4.5)  can  be  rewritten  as 


A»„(ft-)2 
(1  +  H)2*' 


Further,  when  /  varies  over  all  of  W2fc(IRd),  v  varies  over  all  of  £<2(nid),  i.e.,  g  :=  |i/|2  varies  over 
all  non-negative  functions  in  Xi(IRd).  This  means  that  the  ^-approximation  order  requirement  is 
equivalent  to  the  existence  of  c  >  0  such  that 

(4.6)  l»lxc/t(f^|*jyt  <  <*,*lslli.i(B-).  VA  >  0,  ¥5  €  i.OR'1). 


Fixing  h,  the  last  condition  states  that  Xc/h  »  considered  as  a  linear  functional  on  Xj(lRd), 

is  bounded  by  ch2k .  Consequently,  having  {S{4>h)}h  provide  approximation  order  k  is  equivalent 
to  the  existence  of  c  >  0  such  that 

11  A$k(/i-)  „  ^ 

,!(f+ i.1)*i,l~(c/#*)  - ch  * 

The  proof  is  thus  completed,  since  upon  rescaling  the  last  condition  becomes 

(4J)  U(j,+'f.|)Jwc)  <  c. 
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* 


Proof  of  Theorem  1.8.  In  the  case  of  this  theorem,  4>h  =  <t>  for  all  h  >  0.  Using  this  in 
(4.7)  and  letting  h  — ►  0,  we  get  that  (4.7)  is  equivalent  to  |  •  |-fcA0  €  L^C).  6 


Remark.  Note  that  the  cube  C  that  appears  in  the  characterization  of  approximation  orders 
is  entirely  incidental.  Since,  for  every  h ,  is  bounded  by  1,  and  also  (h  +  |  •  |)_fc  is  bounded, 
independently  of  h ,  in  any  complement  of  a  neighborhood  of  the  origin,  the  cube  C  can  be  replaced 
by  any  neighborhood  of  the  origin. 

Another  remark  concerns  the  case  k  =  0  which  will  soon  be  considered  in  the  context  of  density 
orders.  We  have  not  discussed  approximation  order  0  simply  because  of  lack  of  any  mathematical 
interest:  the  requirement  in  this  case  is  vacuous.  This  is  in  agreement  with  Theorem  4.3,  for  the 
boundedness  of  }h  is  also  a  vacuous  condition,  since  each  A^h  is  uniformly  bounded  by  1. 

This  means  that  the  statement  of  Theorem  4.3  is  valid  also  for  k  =  0. 


With  Theorem  4.3  in  hand,  we  turn  our  attention  to  the  characterization  of  density  orders. 
Our  result  concerning  density  orders  is  as  follows. 


Theorem  4.8.  For  {<£&}/,  C  L2(JR.d),  the  sequence  {S(<f>h)}h  provides  density  order  k  if  and  only 
if 


{ 


(a  +  I-D* 


U 


is  bounded  in  L<x,(C),  and 


(4.9) 


lim  h 
fc—o 


-d 


lhaC 


Ak 


(h  +  |  •  IP 


0, 


Vo  >  0. 


Proof.  In  view  of  Theorem  4.3  and  the  definition  of  density  orders,  the  theorem  here  is 
proved  as  soon  as  we  show  that,  under  the  assumption  that  {(/1+*fl|)t  } fa  >s  bounded,  the  condition 


(4.10)  lim  hd'2~kE(U,  5(<M)  =  0,  V/  G  Wk( Dld) 

rl—*0 


is  equivalent  to  (4.9).  For  this  we  can  follow  the  proof  of  Theorem  4.3  up  to  (4.6)  to  conclude  that 
(4.10)  is  equivalent  to  the  condition  that 

(4.n) 

Choosing  g  :=  \aC  *n  (4-11)  and  rescaling,  we  obtain  (4.9),  so  that  the  necessity  of  (4.9)  for 
^-density  order  is  proved. 

To  prove  the  sufficiency,  we  define 


Ah  —  n  xCA(i  +  |.|)2*’ 


h  >  0. 
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We  view  the  X^  as  elements  of  Li(THd)*.  We  want  to  show  that  (4.11)  holds,  namely  that  {A*}* 
converges  weak-*  to  0.  We  know  that  {Xh}h  are  positive,  uniformly  bounded,  and  by  (4.9), 
Xf,(x  )  — ►  0  for  every  a  >  0.  This  latter  condition  implies  that  X hix^)  0  f°r  any  com¬ 
pact  By  linearity,  Xh(g)  tends  to  0  for  each  compactly  supported  simple  function  g.  Since  such 
functions  are  dense  in  /^(ffi/*),  we  obtain  (4.11).  ^ 

Proof  of  Theorem  1.7.  Since  (h  +  |  •  |)"2fc  <  |  •  |"2fc,  (l-8)  '“P^8  *hat 


lim  h~d  f  .  pTtfk  = 

JhC  (A  +  I  *  I) 

which  is  the  case  a  =  1  in  (4.9),  and  implies  the  rest  of  (4.9),  since  here  <j>h  =  <t>  for  all  h,  hence  A* 
does  not  change  with  h.  Thus,  Theorem  4.8  implies  the  sufficiency  of  (1.8). 

On  the  other  hand,  if  S(<f> )  provides  density  order  k ,  then  (4.9)  holds  (with  A,ph  =  A^,  all  h). 
Since  Ij/I-2*  <  c(h  +  |y|)-2*  for  y  €  hC  \  ( hC/2 )  and  some  absolute  constant  c,  we  obtain  from 

(4.9)  (with  a  =  1) 


(4.12) 


A$(y) 

|y|2* 


<  <h)hd 


where  linu— o*(h)  =  0.  Summing  these  estimates  gives 

(4.i3)  /  Sr  i  hd- 


Since  the  right  side  of  (4.13)  is  o(hd),  we  obtain  the  necessity  of  (1.8).  * 

Combining  the  two  last  theorems  with  Theorem  3.12,  we  obtain 

Theorem  4.14.  Let  {<S/J  be  a  sequence  of  shift-invariant  spaces.  For  each  h,  let  4>h  be  the  best 
approximation  from  Sh  to  gm  =  Then,  {Sh}h  provides  approximation  order  k  if  and  only  if 


f 

{(h  +  \\‘\\)k 

is  bounded  in  Loo(C),  and  { Sh}h  are  kth-order  dense  if  and  only  if,  in  addition  to  the  above, 


(4.15) 


lim  h~d  j  jj—r  jrrfiSk  ~  v  a  >  °- 


Proof  of  Theorem  1.9.  This  follows  from  Theorem  1.6,  Theorem  1.7,  and  the  reduction  to 
the  principal  shift-invariant  case  given  by  Theorem  3.12  (with  <f>*h  =  <t>*  =  Ps9*  f°r  ♦ 

5.  The  Strang-Fix  conditions 

As  mentioned  in  the  introduction,  approximation  orders  from  the  scaled  spaces  {Sk}h  were 
characterized  in  [SF]  under  the  assumptions  that  (a)  the  space  Sh  is  obtained  as  the  h-dilate  of 
the  same  principal  shift-invariant  space  <S(0);  (b)  the  generator  <f>  of  S(4>)  is  compactly  supported; 
and  (c)  the  approximation  order  is  realized  in  a  controlled  manner.  The  controlled  approximation 
assumption,  in  turn,  forces  the  condition  &(0)  ^  0. 
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In  order  to  compare  these  conditions  to  the  characterization  of  approximation  orders  for  prin¬ 
cipal  shift-invariant  spaces  that  we  obtain  in  the  present  paper,  we  assume  in  this  section  that  we 
have  in  hand  a  sequence  {S(<f>h)}h  of  principal  shift-invariant  spaces  which  satisfy  one  or  both  of 
the  following  conditions,  in  which  Si  is  some  neighborhood  of  the  origin,  and  rj  and  p  are  positive 
constants: 

(5.1)  3Si,p,h0  s.t.  |^,(x)|  <  ft  a.e.  on  SI,  V  0  <  h  <  h0; 

(5.2)  3  SI,  Tf,  ho  s.t.  rj  <  |^h(x)|  a.e.  on  SI,  V  0  <  h  <  ho. 

Note  that,  in  case  <j>h  does  not  change  with  h  (i.e.,  when  assumption  (a)  above  holds),  and  <f>  is 
continuous  at  the  origin  (e.g.,  <f>  is  compactly  supported,  as  in  assumption  (b)  above),  (5.1)  is 
satisfied  automatically  and  (5.2)  is  reduced  to  the  mere  condition 

(5.3)  SO)  t  0. 

We  recall  (see  the  remark  after  the  proof  of  Theorem  1.6)  that  the  uniform  boundedness  required 
in  Theorem  4.3  for  ^-approximation  order  can  be  checked  in  any  neighborhood  SI  of  the  origin, 
hence  we  can  replace  the  cube  C  in  the  theorem  by  SI.  As  the  next  results  show,  can  often  be 
replaced  by 

(5.4)  M*:=(  Y.  A.(-+0)l’)'/2  =  ([&.?*]  -|W)1/2- 

/3€2irZ*\0 


Lemma  5.5.  If  (5.1)  holds  and  the  sequence  {S(4>h.)}h  provides  approximation  order  k,  then 

(5-6)  t(h  +  \.\)^h<h'° 

is  bounded  in  L^Si')  for  some  0- neighborhood  Si '  and  some  h'0  >  0.  On  the  other  hand,  if  (5.2) 
holds  and  (5.6)  is  bounded  in  Loo(Si')  for  some  0- neighborhood  Si  and  some  h'0,  then  {S(4>h)}h 
provides  approximation  order  k. 


Proof.  If  {S(<t>h.)}h  provides  approximation  order  k,  then,  by  Theorem  4.3, 
{(h  +  |  •  |)-feA^}/i  is  bounded,  say  by  c,  on  Si.  This,  together  with  (5.1),  implies  that 

(5.7)  (5  +  1-  <  c(Mj  +  |5»|!)  <  c(Mj  +  ,■»), 


and  therefore, 

«fc+|.|)-«-c)Mi[<V. 


Thus,  for  sufficiently  small  h  and  some  neighborhood  Si'  C  Si  of  the  origin,  the  leftmost  term  in 
(5.7)  does  not  exceed  2 cp2. 

Conversely,  (5.2)  implies  that,  on  Si, 


|S|2  ^ 

M2h  +  |^|2  ‘  \$h\2 


<  q“2M*. 


Therefore,  by  Theorem  4.3,  the  boundedness  of  (5.6)  implies  that  {S(0j»)}a  provides  approximation 
order  k.  4 

We  now  consider  in  more  detail  necessary  conditions  for  approximation  order  which  follow 

a  » 

from  our  characterization  of  approximation  order.  Since  |^/t(*  +  (3)\  <  Ms  for  all  p  €  2jtZZ  \0,  the 
next  theorem  is  a  direct  consequence  of  the  last  lemma: 
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Theorem  5.8.  If  (5.1)  holds  and  {S{4>h)}h  provides  approximation  orderk,  then,  foraJJO  <  h  <  ho 
and  for  all  (3  €  2w7Zd\0,  and  in  some  O-neighborhood, 

!&.(•  + /?)l<c(*  +  M)fc, 


for  some  c  independent  of  f3  and  h. 

In  case  <f>  does  not  change  with  h,  we  may  let  h  — ►  0  in  the  last  display  and  so  obtain  Theorem 
1.14.  This  shows  that  the  necessity  of  the  Strang-Fix  conditions  (1.11)  for  ^-approximation  order 
holds  in  a  very  general  setting.  This  is  remarkable,  since  this  implication  is  considered  to  be  the 
“harder”  one.  An  analogous  Zoo-result  has  been  obtained  in  [BR2]  by  other  means. 

We  now  consider  in  more  detail  sufficient  conditions  for  approximation  order.  There  is  no 
reason  to  believe  that  (upon  assuming  (5.2))  the  assumptions 

(5.9)  =  0  on  2ir7Ld\0  for  all  M  <  k 

would  suffice  for  approximation  order  k  since  from  Lemma  5.5  we  only  can  deduce  the  following: 
Corollary  5.10.  If  0  <  tj  <  <j>  a.e.  on  some  neighborhood  fl  of  the  origin,  and  if 

(5.11)  ^2  \$(-  +  P)\2  <c\-\2k,  a.e.  on  fl, 

0€2irZd\O 

then  S(4>)  provides  approximation  order  k. 

However,  assumptions  like  (5.9)  can  only  imply  that,  for  each  individual  /?  €  27r2Zd\0, 

hence  will  not  in  general  yield  (5.11).  On  the  other  hand,  there  are  several  results  in  the  literature 
which  show  that,  under  additional  assumptions  on  <t>,  (5.9)  does  imply  that  S(<t>)  provides  approx¬ 
imation  order  k.  For  example,  standard  polynomial  reproduction/quasi-interpolation  arguments 
show  that  if 

(5.12)  M*)|  =  0{\x\~k~d~e),  as  x  oo, 

and  if  <t>{0)  /  0,  then  (5.9)  implies  that  S(4>)  provides  approximation  order  k  (cf.  e.g.,  Proposition 
1.1  and  Corollary  1.2  in  [DJLR]).  Unfortunately,  the  decay  conditions  (5.12)  fail  to  hold  for  many 
functions  <f>  of  interest  (primarily  radial  basis  functions,  and  usually  because  <f>  is  not  smooth 
enough  at  0),  and  in  such  a  case,  the  polynomial  reproduction  argument  either  fails,  or  is  not  easily 
converted  into  approximation  orders.  Circumventing  the  polynomial  reproduction  argument  was 
actually  the  major  objective  of  [BR2].  In  our  context,  Theorem  1.6  leads  to  a  remarkable  result, 
which  allows  (5.12)  to  be  replaced  by  a  much  weaker  condition,  and  which  we  now  describe. 
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For  this  result,  we  need  a  local  version  W2p(ft)  of  the  potential  spaces  If  P  *s  ar> 

integer,  then  this  space  is  simply  the  Sobolev  space  of  all  functions  whose  (weak)  derivatives  up 
to  order  p  (inclusive)  are  in  Li(Q.).  In  this  case,  if  {ftp}pg/  is  a  disjoint  collection  of  open  subsets 
of  IRd,  we  have  SZ/jg/ ll/ll =  ll/li *  As  ‘s  weh-hnown,  there  are  several  equivalent 
extensions  of  the  definition  of  Wf(0.)  to  the  case  of  a  fractional  p  (see,  e.g.,  [A;  ch.7]).  For  fractional 
p,  we  have  the  following  subadditivity  property: 

(5-13)  TjfW'w-'o.*  -  c\\f\\w'(ugeinfi)' 

Pei 

whenever,  say,  {ftp}p  is  a  disjoint  collection  of  cubes;  (cf.  [A;  p.225]).  Our  result  is  as  follows: 

Theorem  5.14.  Assume  that  0  <  i?  <  <f>  a.e.  on  some  cube  ft  centered  at  the  origin.  Let  A  := 
U^€2TZ«i\o(fl  +  0)-  If  <f>  €  Wf(A)  for  some  p  >  k  +  d/2,  and  if  (5.9)  holds,  then  S((f>)  provides 
approximation  order  k. 

The  virtue  of  this  theorem  is  that  we  can  take  ft  to  be  so  small  that  A  does  not  contain 
the  origin.  This  is  important  since  in  many  cases  of  interest  </>  is  smooth  on  IR  \0  but  has  some 
singularity  at  the  origin  (this  happens,  e.g.,  when  <j>  is  obtained  by  the  application  of  a  difference 
operator  to  a  fundamental  solution  of  an  elliptic  equation).  But,  if  4>  satisfies  (5.12),  then  (f>  is 

^  j 

globally  smooth,  since  we  obtain  from  (5.12)  that  </>  G  VF2  (IR  )  for  p  =  k  +  d/2  -f  e/2  as  well  as 
<t>  G  Cfc(IR‘<).  Thus,  Theorem  5.14  and  Theorem  1.14  together  imply  the  following  result. 

Corollary  5.15.  If  <f>  satisfies  (5.12)  and  </>(0)  /  0,  then  S(<fi)  provides  approximation  order  k  if 
and  only  if  (5.9)  holds. 

Proof  of  Theorem  5.14.  It  follows  from(5.9)  that,  for  every  0  G  27rZd\0,  and  with  ftp  :=  ft  +  /3, 

(5.16)  \4>{x  +  /?)|  <  c|i|fc  max||D>||i<)o(ntf),  forx  G  ft. 

m-* 

Since  p  >  k  +  d/2,  the  Sobolev  embedding  theorem  (cf.  [A;  p.217])  implies  that  VF2  (ftp)  is  contin¬ 
uously  embedded  in  the  Sobolev  space  W^ftp).  Thus, 

with  ci  independent  of  (3  (since  all  the  ftp  are  translates  of  each  other).  Substituting  this  into  (5.16) 
we  obtain  that 

\<j>(x  -f  0) |  <  C2|x|fc||^||iv^(ofl)»  x  €  ft,/3  G  2ir7Ld\0. 

Squaring  the  last  inequality  and  summing  over  0  G  2v7Zd\0,  we  obtain,  in  view  of  (5.13),  that 

P&irTL^O 

Lemma  5.5  now  supplies  the  conclusion  that  S(<f>)  provides  approximation  order  k.  4 


18 


In  applications,  it  might  be  convenient  to  take  p  to  be  the  least  integer  that  satisfies  p  >  k+d/ 2. 

For  this  case,  Theorem  5.14  reduces  to  Theorem  1.15. 

References 

[A]  R.A.  Adams,  Sobolev  Spaces ,  Academic  Press,  New  York,  1975. 

[Bl]  C.  de  Boor,  The  polynomials  in  the  linear  span  of  integer  translates  of  a  compactly  supported 
function,  Constructive  Approximation  3  (1987),  199-208. 

[B2]  C.  de  Boor,  Quasiinterpolants  and  approximation  power  of  multivariate  splines.  Computation 
of  curves  and  surfaces,  M.  Gasca  and  C.  A.  Micchelli  eds.,  Dordrecht,  Netherlands:  Kluwer 
Academic  Publishers,  (1990),  313-345. 

[Bui]  M.  D.  Buhmann,  Multivariate  interpolation  with  radial  basis  functions,  Constructive  Approx¬ 
imation  0  (1990),  225-256. 

[Bu2]  M.D.  Buhmann,  On  quasi-interpolation  with  radial  basis  functions,  ms.  (1991). 

[BD]  C.  de  Boor  and  R.  DeVore,  Approximation  by  smooth  multivariate  splines,  Transactions  of 
Amer.  Math.  Soc.  276  (1983),  775-788. 

[BuD]  M.D.  Buhmann  and  N.  Dyn,  Error  estimates  for  multiquadric  interpolation,  in  Curves  and 
Surfaces,  P.J.  Laurent,  A.  Le  Mehaute  and  L.L.  Schumaker  (ed.),  Academic  Press,  New  York, 
1990,  1-4. 

[BH1]  C.  de  Boor  and  K.  Hollig,  B-splines  from  parallelepipeds,  J.  d’Anal.  Math.  42  (1982/3), 
99-115. 

[BH2]  C.  de  Boor  and  K.  Hollig,  Approximation  order  from  bivariate  C^-cubics:  a  counterexample, 
Proc.  Amer.  Math.  Soc.,  87  (1983),  649-655. 

[BJ]  C.  de  Boor  and  R.Q.  Jia,  Controlled  approximation  and  a  characterization  of  the  local  approx¬ 
imation  order.  Proc.  Amer.  Math.  Soc.,  95  (1985)  547-553. 

[BR1]  C.  de  Boor  and  A.  Ron,  The  exponentials  in  the  span  of  the  integer  translates  of  a  compactly 
supported  function:  approximation  orders  and  quasi-interpolation,  J.  London  Math.  Soc.,  to 
appear. 

[BR2]  C.  de  Boor  and  A.  Ron,  Fourier  analysis  of  approximation  orders  from  principal  shift-invariant 
spaces,  ms. 

[C]  C.  K.  Chui,  Multivariate  Splines,  CBMS-NSF  Regional  Conference  Series  in  Applied  Math., 
#54,  SIAM,  Philadelphia,  1988.  (1935),  pp.vii+152. 

[CL]  E.W.  Cheney  and  W.A.  Light,  Quasi-interpolation  with  basis  functions  having  non-compact 
support,  Constructive  Approx.,  to  appear. 

[DJLR]  N.  Dyn,  I.R.H.  Jackson,  D.  Levin,  A.  Ron,  On  multivariate  approximation  by  the  integer 
translates  of  a  basis  function,  Israel  J.  Math.,  to  appear. 

[DM1]  W.  Dahmen  and  C.  A.  Micchelli,  Translates  of  multivariate  splines,  Linear  Algebra  and  Appl. 
52/3  (1983),  217-234. 

[DM2]  W.  Dahmen  and  C.A.  Micchelli,  On  the  approximation  order  from  certain  multivariate  spline 
spaces,  J.  Austral.  Math.  Soc.  Ser  B  26  (1984),  233-246. 


19 


[DR]  N.  Dyn  and  A.  Ron,  Local  approximation  by  certain  spaces  of  multivariate  exponential- 
polynomials,  approximation  order  of  exponential  box  splines  and  related  interpolation  prob¬ 
lems,  Transactions  of  Amer.  Math.  Soc.  319  (1990),  381-404. 

[H]  H.  Helson,  Lectures  on  Invariant  Subspaces ,  Academic  Press,  New  York,  1964. 

[HL]  E.  J.  Halton  and  W.  A.  Light,  On  local  and  controlled  approximation  order,  preprint  1990. 

[J 1]  R.Q.  Jia,  Approximation  order  from  certain  spaces  of  smooth  bivariate  splines  on  a  three- 
direction  mesh,  Trans.  Amer.  Math.  Soc.  295  (1986),  199-212. 

[J2]  R.Q.  Jia,  A  counterexample  to  a  result  concerning  controlled  approximation,  Proc.  Amer. 
Math.  Soc.  97  (1986),  647-654. 

[JL]  R.Q.  Jia  and  J.J.  Lei,  Approximation  by  multiinteger  translates  of  functions  having  non¬ 
compact  support,  J.  Approximation  Theory,  to  appear. 

[Ja]  I.R.H.  Jackson,  An  order  of  convergence  for  some  radial  basis  functions,  IMA  J.  Numer.  Anal. 
9  (1989),  567-587, 

[L]  W.A.  Light,  Recent  developments  in  the  Strang-Fix  theory  for  approximation  orders,  Curves 
and  Surfaces,  P.J.  Laurent,  A.  le  Mehaute,  and  L.L.  Schumaker  eds.,  Academic  Press,  New 
York,  1990,  xxx-xxx. 

[LJ]  J.J.  Lei  and  R.Q.  Jia,  Approximation  by  piecewise  exponentials,  SIAM  J.  Math.  Anal.,  to 
appear 

[M]  W.  R.  Madych,  Error  estimates  for  interpolation  by  generalized  splines,  preprint. 

[MN1]  W.R.  Madych  and  S.A.  Nelson,  Polyharmonic  cardinal  splines,  I,  J.  Approximation  Theory, 
40  (1990),  141-156. 

[MN2]  W.R.  Madych  and  S.A.  Nelson,  Multivariate  interpolation  and  conditionally  positive  functions 
II,  Math.  Comp.  54  (1990),  211-230. 

[P]  M.  J.  D.  Powell,  The  theory  of  radial  basis  function  approximation  in  1990,  DAMTP  rep. 
NAll,  University  of  Cambridge. 

[R]  A.  Ron,  A  characterization  of  the  approximation  order  of  multivariate  spline  spaces,  Studia 
Math.,  to  appear. 

[Ra]  C.  Rabut,  Polyharmonic  cardinal  B-splines,  Parts  A  and  B,  preprint  (1989). 

[RS]  A.  Ron  and  N.  Sivakumar,  The  approximation  order  of  box  spline  spaces,  Proc.  Amer.  Math. 
Soc.,  to  appear. 

[S]  I.  J.  Schoenberg,  Contributions  to  the  problem  of  approximation  of  equidistant  data  by  analytic 
functions,  Parts  A  &  B,  Quarterly  Appl.  Math.  IV  (1946),  45-99, 112-141. 

[SF]  G.  Strang  and  G.  Fix,  A  Fourier  analysis  of  the  finite  element  variational  method.  C.I.M.E.  II 
Ciclo  1971,  in  Constructive  Aspects  of  Functional  Analysis,  G.  Geymonat  ed.,  1973,  793-840. 


20 


